Dynamical Regimes Induced by Spin Transfer in Magnetic Nanopillars. 
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We demonstrate the predicted out-of-plane precession induced by spin transfer in magnetic nanos- 
tructures with in-plane magnetic field. We show that other magnetic excitations have a significant 
effect on the stability of the out-of plane precession, making it extremely sensitive to the orienta- 
tion of the applied magnetic field. The data are supported with micromagnetic simulations. Our 
results elucidate the relation between the excitation spectrum and the specific dynamical behaviors 
of nanoscale magnets. 
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Spin transfer torque (ST) [l[ exerted on nanomag- 
nets by spin-polarized current / can induce dynamical 
states not accessible by any other techniques, providing 
a unique opportunity to test our understanding of nano- 
magnetism [2j. At / just above the excitation threshold, 
ST acting on the magnetic moment m\ of a nanomagnet 
Fi causes precession on an elliptically shaped orbit, de- 
termined by a combination of the magnetic field H and 
the anisotropy of Fi. The nanomanget is usually a thin 
film whose anisotropy is dominated by the demagnetizing 
field, with H commonly in the film plane. The precession 
amplitude grows with /, resulting in an orbit known as 
the clamshell mode (left inset in Fig. [1} . Good agreement 
between calculations and experiments has been achieved 
for precessional dynamics in this regime [H, Q| ■ 

However, fast operation of magnetic devices driven by 
ST will likely be performed in the less explored high- 
current regime. According to simulations 0, H, @], the 
extreme points of clamshell eventually merge, resulting 
in a crossover to the out-of-plane (OP) precession mode 
consisting of either the lower or the upper half of the 
clamshell. Only a broad spectral feature indicative of 
this mode has been seen suggesting that micromag- 
netic simulations may not be adequate for the highly ex- 
cited dynamical states of nanomagnets, or that current- 
induced effects may not be fully described by the estab- 
lished ST mechanisms. 

ST also affects the magnetic layer F2 used to polarize 
the current, which can decrease the dynamical coherence 
and suppress the OP mode p}. Here, we report observa- 
tion of the OP precession in devices where this effect was 
minimized by using an extended F2. We identified and 
analyzed the effects of varied direction and magnitude 
of H . The micromagnetic simulations support our inter- 
pretation of the data, and provide insight into the mi- 
croscopic origins of the observed behaviors. Some of the 
results could not be reproduced by simulations, suggest- 
ing that the understanding of current-induced dynamics 
in nanomagnets is still incomplete. 

Multilayers Cu(50)Py(20)Cu(5)Py(5)Au(20), where 
Py=NigoFe2o and thicknesses are in nm, were deposited 
on oxidized silicon at room temperature (RT) by mag- 



netron sputtering at base pressure of 5 x 1CP 9 Torr, in 
5 mTorr of purified Ar. Fi=Py(5) and about 5 nm- 
thick part of F2=Py(20) were patterned by Ar ion 
milling through an evaporated Al mask with dimensions 
of fOO nm x 50 nm, followed by deposition of 30 nm of 
undoped Si without breaking the vacuum. This proce- 
dure avoids oxidation of the magnetic layers, which can 
affect the magnetic dynamics 8] . The mask was removed 
by a combination of ion milling with Ar beam nearly par- 
allel to the sample surface, and etching in a weak solution 
of HF in water, followed by sputtering of a 200 nm thick 
Cu top contact. We discuss data for one of three devices 
that exhibited similar behaviors. 

All measurements were performed at RT. The sample 
was contacted by coaxial microwave probes, which were 
connected through a bias tee to a current source, a lock- 
in amplifier, and a spectrum analyzer through a broad- 
band amplifier. To enable the detection of the preces- 
sional states by electronic spectroscopy, H was rotated 
in the sample plane by angle <p = 40° with respect to 
the nanopillar easy axis, unless specified otherwise. Pos- 
itive I flowed upwards. The device was characterized by 
magnetoresistive (MR) measurements of its response to 
H and /, yielding the parameters essential for modeling, 
such as the MR of 0.2f fi, the dipolar coupling field of 
200 Oe caused by the partial patterning of F2, and the co- 
ercivity of 175 Oe. The latter was consistent with Stoner- 
Wohlfarth approximation, indicating uniform magnetic 
reversal. The measured microwave signals were adjusted 
for the frequency-dependent gain of the amplifier and 
losses in the cables and probes, determined with a cali- 
brated microwave generator and a power meter. 

Fig. [T] shows the dependence of the measured power 
spectral density (PSD) on H. The spectra at H — 550 Oe 
exhibit three harmonically related peaks at / < 5 mA, 
caused by the clamshell precession. The expected fun- 
damental frequency of the OP mode is close to the 
frequency of the clamshell's second harmonic near the 
crossover, since its trajectory is half that of clamshell. 
The 550 Oe data in Fig. Q] exhibit a spectral feature at 
I > 5 mA consistent with this relationship between the 
two modes. Below, we present measurements and micro- 
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FIG. 1: PSD vs. frequency / and /, at the labeled values of 
H. The logarithmic data scale brings out a break at 2.93 GHz 
of about 50 fW /MHz due to the spectrum analyzer crossover. 
Inset: clamshell trajectory of the magnetic moment, with H 
shown. 



magnetic simulations that confirm our interpretation and 
elucidate other, more complex dynamical behaviors. 

At H > 550 Oe, the OP peak rapidly shifts to higher 
frequency. Simultaneously, the high-current part of the 
peak splits from the lower part, and eventually merges 
with the clamshell peaks (H = 650 Oe and 700 Oe 
data in Fig. [JJ. One can estimate the precession am- 
plitude based on the total microwave power under the 
peaks divided by I 2 (see also the description of sim- 
ulations). The largest possible emitted power for hy- 
pothetical oscillations between the parallel (P) and an- 
tiparallel (AP) configurations of the magnetic layers is 
37.7 pW/mA 2 for our sample. The power generated 
by clamshell precession at H = 700 Oe, I — 3.5 mA 
is 11.2 pW/mA 2 , and by the OP mode at the same H 
and / = 10 mA is 9.9 pW/mA 2 . Both values correspond 
to the in-plane precession angle exceeding 90°, providing 
a strong evidence for our interpretation of the spectral 
peaks as large-angle precessional modes. If these spec- 
tral features were induced by inhomogeneous dynamics 



rather than precession, they would result in significantly 
smaller microwave power emission. The first clamshell 
harmonic exhibits the smallest FWHM of 30 MHz at 
H = 550 Oe, / = 3.5 mA. The peaks decrease in am- 
plitude and broaden with increasing H. These behaviors 
suggest increasingly inhomogeneous dynamics, resulting 
in decoherence of precession. In contrast, the intensity of 
the OP peak at I = 10 mA increases from 0.86 pW/MHz 
at H = 550 Oe to 2.5 pW/MHz at 1 kOe, while the 
FWHM remains approximately constant at 140±15 MHz. 

The spectra are asymmetric with respect to the direc- 
tion of H , as illustrated by the difference between the left 
and the right panels in Fig. [TJ The clamshell peaks are 
consistently broader at H < than at H > 0. At large /, 
they are replaced by incoherent noise rather than the OP 
mode. The direction of H for which the OP mode was 
observed varied among the samples, indicating extrin- 
sic origin of asymmetry. At large H, the data became 
similar for positive and negative H. In particular, the 
H = ±1 kOe data exhibit a sharp OP peak at I > 9 mA. 

To gain insight into the origin of the behaviors shown in 
Fig. [TJ as a well as other results discussed below, we per- 
formed micromagnetic simulations with OOMMF open 
code software [9j. The simulations included the current- 
induced ST and Oersted field effects, but neglected ther- 
mal fluctuations. The dipolar field was accounted for 
by subtracting 200 Oe from the in-plane component of 
H. The cell size was 4x4x5 nm 3 . Reducing the size 
to 2 x 2 x 5 nm 3 did not significantly affect the results. 
The parameter values available mostly from MR mea- 
surements (Io| were used: current polarization p = 0.7, 
Py exchange stiffness A — 1.3 x 10 -6 erg/cm, Gilbert 
damping a = 0.03, and the ratio A = 1.3 of the ST mag- 
nitudes near the AP and the P states. Saturation mag- 
netization M = 750 emu/cm 3 of Py was determined by 
magnetometry of a Py(5) film prepared under the same 
conditions as the nanopillar. 

The spectra were calculated from the simulated time- 
dependent magnetization distribution of the nanopillar. 
The calculated time-dependent resistance was R(t) — 
R + AR(t), where R = (R P + R A p) /I, and AR(t) = 
(Rp — Rap){si -s 2 )/2. Here, s 1 (t),s 2 (t) are the local 
normalized magnetizations of Fi and F2, and () de- 
notes averaging over the simulation grid. The ac volt- 
age on the input of the amplifier V(t) = i^^'^qq wa $ 
obtained by assuming that a constant current / is dis- 
tributed between a 50 f2 load and the sample resistance 
R(t). Fast Fourier transform (FFT) of V(t) over a pe- 
riod of T — 16.4 ns with a 1 ps step was performed 
after relaxation for 10 ns. The power spectral density 
was determined by PSD(f) = 2V 2 (f)/(50flAf), where 
A/ = 1/T, and a factor of 2 accounts for the negative-/ 
contribution to the FFT. 

Our calculations showed that a significant asymmetry 
of spectra with respect to the direction of H can be in- 
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FIG. 2: (a) Calculated PSD for # = 650 Oe. (b) same 
as (a), for H = —650 Oe. c, instantaneous magnetization 
distribution for H — —650 Oe, 1 — 6 mA, captured in the 
I-II-III sequence with a 100 ps interval, showing the nucle- 
ation, propagation, and annihilation of a vortex; the mo- 
tion of its core is marked with yellow dots. The intensity 
of blue(red) reflects the out-of-plane component of the mag- 
netization above(below) the plane. Arrow shows the direction 
of H. 



duced by the simultaneous effects of the current-induced 
Oersted field and a modest asymmetry of the sample ge- 
ometry. Fig. [5] shows results for a nanopillar approxi- 
mated by two semi-ellipses with minor semi-axes of 34 nm 
and 18 nm, and a major axis of 104 nm. Coherent OP 
precession was obtained at H > 0, but was suppressed at 
H < by vortices and finite- wavelength spin- waves 11 1 . 
The vortices usually nucleated at the right upper edge 
of the nanopillar, and annihilated at the lower left edge 
(Fig. • The chirality of the vortices coincided with the 
direction of the Oersted field's rotation, indicating that 
the asymmetry of spectra is caused by the suppression or 
enhancement of vortex nucleation due to the interplay of 
sample shape asymmetry and the effect of Oersted field. 

Simulations could not reproduce several features of the 
data for any reasonable variations of nanopillar shape, 
distribution of current and its polarization, or Py stiff- 
ness. Firstly, the simulated OP peak did not exhibit 
the rapid shift and splitting with increasing H seen in 
data. Secondly, despite a significant asymmetry of the 
calculated spectra, they did not reproduce the region 
at / > 5 mA where sharp OP peak was present for 
H = 500 Oe, but no dynamical features appeared at 
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FIG. 3: Frequency of the OP peak vs I for (f> = 48° (squares), 
4> = 40° (open circles) , (j> = 24° (triangles) . (a) Measurements 
performed at H = 550 Oe, (b) Calculations performed at 
H = 650 Oe, with the same notations. 



H = —550 Oe. The simulations also indicated that 
the OP mode should exhibit multiple spectral harmonics, 
while only one or two harmonics could be seen in data, re- 
gardless of the large amplitude of precession established 
from the analysis of peak intensity. These features may 
be caused by additional effects of spin-polarized current 
neglected by the model, or by the dynamical states not 
described by micromagnetic simulations. 

The differences between data and simulations also open 
the possibility that the high-current spectral feature in 
Fig.[T]is caused by dynamics different from the OP mode. 
Current-induced excitation of the polarizing layer F2 can 
lead to microwave peaks, appearing above the onset cur- 
rent determined by the ratio of the volumes of F2 and 
Fi 0]. However, the effective volume of the extended 
layer F2 in our samples far exceeds that of Fx, and thus 
cannot explain the onset current that at 550 Oe is only 
2.6 times larger than the onset of the clamshell preces- 
sion. The nanopillar shape imperfections can also re- 
sult in precession around a configuration intermediate 
between the AP and P states. However, this intermediate 
state would quickly become unstable at increased H, con- 
trary to the high-field data in Fig. Q] Such a state would 
also likely appear as an intermediate-resistance step not 
seen in our dc measurements of MR. 

Both measurements and simulations showed that the 
OP precession is extremely sensitive to the orientation 
of H , when the latter was rotated in the film plane or 
tilted out of plane. Rotating H in the plane changed the 
dependence of the OP mode frequency on /, Fig. (3ja). 
At 4> = 48°, the peak exhibited a blue shift up to 6.5 mA, 
above which it gradually red shifted. At smaller values 
of </>, the peak broadened, decreased in intensity, and red 
shifted. The correlation between the width of the OP 
peak and the dependence of its frequency on / is con- 
sistent with published simulations [2, 0]- Namely, blue 
shift is always predicted in the macrospin approximation, 
which is more applicable to narrow coherent peaks. In 
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FIG. 4: PSD for \H\ = 900 Oe tilted by 45° with respect to the 
sample plane, at <f> = 40°. Insets are side view schematics for 
the orientation of H with respect to the pillars, where positive 
in-plane direction of H is to the right, (a)-(d) Measurements, 
(e)-(h) Simulations. 



contrast, the peaks can red shift in micromagnetic simu- 
lations of the more inhomogeneous dynamics associated 
with broader spectral peaks. The red shift originates 
from the decrease of the total magnetic moment of the 
nanopillar caused by the inhomogeneity. This interpreta- 
tion was supported by our simulations (Fig.[3Jb)), where 
decreasing <fr resulted in increasingly inhomogeneous OP 
dynamics. The OP mode red shifted with / at small <fi. 
At larger (f>, it blue shifted at small / and red shifted at 
large /, in excellent overall agreement with the data. A 
somewhat larger H was used in simulations to destabi- 
lize the static current-induced AP state at large I (see 
discussion of Fig. [J). 

For H tilted with respect to the sample plane, the dy- 
namics showed dependence not only on the sign of H , but 
also on its tilting direction (Fig. 2]). The component of H 
perpendicular to the film plane was significantly smaller 
than the demagnetizing field of 9.4 kOe, resulting in only 



5° tilting of mi when it was static. The OP mode ap- 
peared for both directions of H above the film plane, and 
was suppressed for H below the film plane. The different 
response of dynamics to the opposite directions of tilt- 
ing indicates a significant asymmetry of this state with 
respect to the film plane. This asymmetry characteris- 
tic of the OP dynamics was reproduced by the simula- 
tions, as shown in panels (e)-(h). Simulations for H tilted 
above the sample plane showed predominantly coherent 
OP precession below the plane, in excellent agreement 
with robust spectral features seen at large I in the data 
of Figs.|4jb,c). We note that simulations with in-plane H 
also showed OP precession below the plane, determined 
by the relative orientations of H and the nanopillar easy 
axis. In contrast, simulations with H tilted below the 
sample plane showed OP precession above the sample 
plane suppressed by dynamical inhomogeneities. These 
results are consistent with the data. 

In summary, coherent current-induced dynamics in 
magnetic nanopillars was achieved by employing an ex- 
tended polarizing layer. A high-current spectral peak 
was identified as the out-of-plane precession, whose de- 
pendence on the orientation of H was in excellent agree- 
ment with micromagnetic simulations. The dynamics ex- 
hibited a significant asymmetry with respect to the di- 
rection of H , attributed to a combination of nanopillar 
shape asymmetry and Oersted field. 
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